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Abstract. Context. We investigate mid-infrared and X-ray properties of the dusty torus in unification scenarios for active 
galactic nuclei. 

Aims. We use the relation between mid IR and hard X-ray luminosities to constrain AGN unification scenarios. 
Methods. With VISIR at the VLT, we have obtained the currently highest angular resolution ((y.'35 FWHM) narrow-band mid 
infrared images of the nuclei of 8 nearby Seyfert galaxies. Combining these observations with X-ray data from the literature we 
study the correlation between their mid IR and hard X-ray luminosities. 

Results. We find that the rest frame 12.3 jim (Lmir) an d 2-10 keV (L x ) luminosities are correlated at a highly significant level. 
The best fit power-law to our data is log L M i R oc ( 1 .60 ± 0.22) log L x , showing a much smaller dispersion than earlier studies. 
Conclusions. The similarity in the log Lmir / log Lx ratio between Sy Is and Sy2s even using high angular resolution MIR data 
implies that the similarity is intrinsic to AGN and not caused by contamination from extra-nuclear emission. This supports 
clumpy torus models. The exponent of the correlation constrains the inner geometry of the torus. 
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1. Introduction 

The unification model for active galactic nuclei (AGN) inter- 
prets the different appearance of Seyfert 1 (Syl) and Seyfert 
2 (Sy2) galaxies uniquely as the result of an orientation effect 
jAntonuccill 19931 lBarthelll994t lUrrv & Padovanil H996). The 
central engine is considered to be surrounded by an optically 
and geometrically thick molecular torus. Associated with this 
torus are large masses of dust that supposedly reprocess the X- 
ray and UV emission from the accretion disk and re -emit it in 
the mid infrared (MIR) regime dPier & Krolikll993l) . 

It is thus very attractive to search for correlations between 
IR continuum and hard X-ray emission in order to test the 
unification scenarios for AGN. A tight correlation between 
the 10.5 fim continuum and the absorption-corrected 2-10 keV 
luminosi ties for 8 nearby Seyfert galaxies was reported by 
IKrabbe et all fcOOll) usin g 172 resolution MIR data. More re- 
cently |LutzetalJ (|2004) found a correlation between the rest 
frame 6 fim luminosity and the absorption-corrected hard X- 
ray luminosity for a sample of 71 AGN. This sample was com- 
prised of objects for which 24" angular resolution ISOPHOT 
spectra and hard X-ray observations were available; in partic- 
ular it does not contain Compton-thick objects. However, the 
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authors reported two problems that the unification model faces 
when compared with their observations. (I) The scatter of the 
relation is about a n order of ma gnitude larger than expected 
from the results of Krabbe et al. (2001J). (II) There is no sig- 
nificant difference between early type (Sy types 1 to 1.5) and 
late type (Sy types 1.8 to 2) objects in the average ratio of mid- 
infrared to hard X-ray emission, as would be expected from an 
optically and geometrically thick torus dominating the mid-IR 
AGN continuum. 

Here we study this relation through res t wavele ngth 
12.3 yu m p hotometry of AGN , improving upon IKrabbe et alJ 
( 200U) and lLutzetaP ( 12004 by a factor of 3 and 80, respec- 
tively, in terms of angular resolution. The wavelength was cho- 
sen as, among the bands accessible from the ground, it best 
traces hot and warm dust. 

Throughout this pa per we assume H o = 72 km/s, Q.\ = 
0.73 and O m = 0.27 JSpergel et al.l2003l) . 

2. Observations and data analysis 

Between March and July 2005 we observed 8 nearby Seyfert 
nuclei (see Tab le [0, select ed from the samp le used by 
iLutz et alJ ll2004l) . with VISIR llLagage et alJ2004l) at the VLT. 
The individual objects were selected as to avoid a strong bias 
toward higher luminosities with higher redshift. All observa- 
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tions were carried out under good photometric and optical see- 
ing conditions (< 078) and with airmass below 1.3. For all ob- 
servations we used the parallel chop/nod mode of VISIR with a 
chopping throw of 8", thus receiving all three beams on the de- 
tector. We observed each object in three different narrow-band 
filters in order to obtain the rough shape of its MIR spectral en- 
ergy distribution (SED). These filters were chosen as to avoid 
emission and absorption lines, in particular the 12.8 jum Nell 
line and the broad 9.7 /jm silicate feature whic h in cases of very 
strong absorption can extend to A * 12.5 fim dPorschner et alJ 
119781) . NGC 7314 was only observed in two different niters. 
Total integration times ranged between 1.2 and 9.2 min per fil- 
ter and per object. 

We reduced science and standard star frames using the 
pipeline written by Eric Pantin (private communication). To 
eliminate glitches, the pipeline applies a bad pixel mask and re- 
moves detector stripes. Subsequently we removed background 
variations using a 2 dimensional polynomial fit. Here we 
treated each nodding cycle separately as the background pat- 
tern sometimes changes between two consecutive cycles. The 
count rate for one full exposure was calculated as the mean of 
all 3 beams from all nodding cycles of this exposure. As an er- 
ror estimate we use the standard deviation of these. In order to 
minimise the effect of residual sky background we chose rela- 
tively small apertures (=s 10 pixels = 1727) for the photometry 
and corrected the obtained count rates using the radial profiles 
of standard stars. Finally we calibrated our photometry using 
the same standard stars. These were observed within 1.5 hrs 
of the science exposures and with a maximum airmass differ- 
ence of 0.15. Our observational results are shown in Table|2] 
Comparing the FWHM of standard star and science exposures 
yields no indication that our sources are extended. Thus our 
observations indeed probe the physical scales given in Tabled 



3. Results 

Figure ^ shows all photometric data points taken during 
our observing programme. For comparison Fig. [2] shows our 
measurements to gether with ISOPHOT spectra taken from 
Rigppou lou et alJ d 19991) for these three sources common to 
both samples. Among our data the SED of PKS 2048-57 is 
the only one displaying a strong 9.7 jjm silicate feature which, 
however, does not significantly affect our measurement of the 
continuum flux in the Nellrefl filter. This is important as it 
validates our result on the MIR - hard X-ray correlation (see 
below). 

One of the objects in our sample - NGC 4579 - has peculiar 
MIR propert i es. Co mparing the ISOPHOT 6 /im flux reported 
in lLutz et alJ d2004 of S 6fim = 96 mJy to our VISIR flux levels 
(see Table|2ji leads to the conclusion that either this object has a 
peculiar MIR SED that does not rise toward longer wavelengths 
as expected for typical AGN or that the bulk of IR emission 
(> 70%) does not originate in the central point source. In both 
cases it is highly questionable whether the lat ter is dominated 
by a molecular torus. Moreover lContinil d2004h reports the MIR 
emission of NGC 4579 to be dominated by thermal emission 
of shock-heated ISM in the vicinity of the nuclear region rather 



Table 1. Basic parameters of our sample of galaxies. 
Heliocentric redshifts z were taken from th e NED, Seyfert 
types are according to IVeron-Cettv & Veronl J2000l) . The ab- 
sorption corrected hard X-ray luminosities we compiled from 
the literature. Where multiple observations were available the 
average was taken. In these cases the reported uncertainties cor- 
respond to 1/2 of the peak-to-peak variabilities. References are 
^Turner & Pounds! dl989l) . ^Gondoin et alJ d200ll). ^Lutz et all 
d2004. ^Reeves et al.1 J2004, ferashima etaD d2002h 
IHo et alJ (Eool, ^Reynolds et alJ (Eool ^Steenbrugge et alJ 
d2003[) . ICappietal.1 d2005l) ^Risalitil d2Q02l). lLandietalJ 
(boOlL "iKraemer et aH d2004 . " iRisaliti et al.1 h00(ti . 
° ISmith & DoneMl996l) . In the case of PG 2130+099 only one 
measurement could be found. Here we assume a variability 
of a factor of 2. The last column contains the physical scales 
resolved by VISIR, computed for an angular resolution of 
0".35, which is a typical value for our observations (see table 
|2}. 'For a discussion of the classification of NGC 4579 see 
section |3] 
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than by reprocessed X-ray and UV emission from the central 
engine itself. These peculiarities plus the fact that NGC 4579 
has also been cla ssified as a lo w-ionis ation nuclear emission re- 
gion (LINER) bv lKeelld 19831) and lFilippenko & Sargendd 19851) 
led us to exclude it from the MIR - hard X-ray correlation in- 
vestigated here. NGC 7314, the other low luminosity object in 
our sa mple, harbours a regular Seyfert nucleus dLumsden et alJ 
120041) . 

For all objects in our sample we calculated the 12.3 fim rest 
frame luminosity by using the filter with the central wavelength 
closest to 12.3 //m and assuming a flat SED in this spectral re- 
gion. Then we compiled absorption-corrected 2-10 keV mea- 
surements for our sample from the literature (see Tabled for 
details). When several such measurements were available we 
use the mean luminosity for the MIR - hard X-ray correlation 
and take half of the peak-to-peak-variability as a conservative 
estimate of the uncertainty in the hard X-ray luminosity. Only 
one such measurement could be found for PG 2130+099 for 
which we therefore assumed a variability of a factor of 2. As 
we derive the luminosity distances of all targets solely from 
their redshifts we assume an uncertainty of 3 Mpc due to possi- 
ble deviations from the Hubble flow. This is applied to the error 
bars of both X-ray and MIR luminosities. 

The resulting MIR - hard X-ray correlation we obtain for 
our VISIR sample is shown in Fig. [3] Pearson's coefficient is 
r = 0.98 meaning that the correlation is highly significant. The 
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Table 2. Basic observational parameters for all data used: Object name, observing date (in 2005), filter, flux and full width half 
maximum for standard star (STD) and science observations. As error estimates we use the standard deviation of all contributing 
beams (see section|5]for details). The FWHM are the average of all beams of one exposure. All science and corresponding STD 
exposures were obtained within a maximum airmass difference of 0.15. "The central wavelengths for the individual filters are 
9.82 /im for SIVrefl, 10.49 /mi for SIV, 11.25 pan for PAH, 11.88/mi for PAH2ref2, 12.27 pan for Nellrefl and 12.81 pirn for 
Nell. Their transmission curves are shown in the VISIR User's manual at http://www.eso.org/instruments/visir/ 
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Fig. 1. The MIR SEDs of our sample of galaxies. Early- and late-type Seyferts are shown in the left- and right-hand panel 
respectively. Horizontal error bars show the band width of the used filters. The flux densities of NGC 7314 have been offset by 



+ 1-10 lu erg s 1 cm 2 for clarity. 
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Fig. 2. Comparison of ISOPHOT-S spectra (dashed histograms) taken from Rigopo ulou et all i 19991) and our VISIR photometry 
(filled circles). Horizontal error bars show the band width of the filters used. 



best fit power-law to our data is logL M i R = (-25.4 + 9.5) + 
(1.60 + 0.22) ■ logLx with measurement errors considered in 
the fit. Including NGC 4579 into the fit yields an exponent of 
~ 1.28 rather than 1.6 and r = 0.95. 



4. Discussion and conclusions 

We have presented new high-resolution MIR observations of 8 
nearby Seyferts, obtained with VISIR at the VLT. We found 
that their rest-frame 12.3jum and absorption-corrected 2-10 
keV luminosities are strongly correlated following Lmir L^ 6 . 
Early- and late-type Seyfert galaxies follow the same relation. 
At the low-luminosity end the correlation is not well defined 
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log A,L^ at 12.3 io.m [erg/s] 



Fig. 3. Correlation of MIR to absorption-corrected hard X-ray 
luminosities for our VISIR sample. Filled diamonds are early- 
type Seyferts (Type Sy 1 .5 or smaller) while filled squares are 
late-type Seyferts. The power law displayed is the best fit ob- 
tained when omitting NGC 4579 (see section|3]for details). 



as one of the two object's MIR radiation is very probably not 
dominate d by a dusty t o rus. 

Both iKrabbe et alJ fcOOlll and Lutz et alJ d2004 found a 
correlation close to Lmir k L^. The exponent different to the 
one we find is most probably due to the difference in angu- 
lar resolution (24" for ISOPHOT, 1V2 for Krabbe et al.'s sam- 
ple and ~ 0735 for VISIR) and to the choice of the rest frame 
wavelength as the intervening medium is expected to be more 
optically thin at 12.3 //m than at shorter wavelengths. In par- 
ticular, using the corresponding ISOPHOT 6 /mi data for our 
sample also yields an exponent close to 1 . In additio n we find a 
dispersion simila r to the one in Krabbe et alJ the larger 

one reported by iLutz et all d2004h can be attributed to their 
relatively poor angular resolution and the fact that the emis- 
sion at 6 is slightly mo re dependent on inclination effects 
(e.g. see Honi g et alJl2006l) . On the other hand we confirm the 
absence of a significant difference in the Lmir / Lx ratio be- 
tween early and late-type objects found by both other studies. 
At our high angular resolution this is clearly in disagreement 
with smooth torus models, in which a flux difference of one 
order of magnitude is expected for typ e 1 and type 2 AGN 
of the same bolometric luminosity (e.g. lPier & Krolikl ll993). 
This cannot be accounted for by dust clumps situated in the 
ionisation cone (Efstat hiou et alj["l995l) because their contri- 
buti on to the n uclear MIR flux should not ex ceed 50% (see 
Gallia no et al .120051) . Clumpy torus models (e.g. lNenkova et all 
2002) do not predict such an offset and seem to be compatible 
with our results. 

iHonig et all {2006 ) recently used 3D radiative transfer 
modelling to demonstrate that clumpy tori can appear as op- 
tically thin in the MIR with most of the radiation originating 
in the innermost part of the torus. In their model individual 
clouds are optically thick but their filling factor is very small. 
Consequently lines-of-sight toward the innermost clouds are 
usually not blocked by other clouds. This implies that Lmir 
is proportional to the surface area A of the inner edge of the 
torus and the covering factor of clouds projected on A. Our 



finding that Syls and Sy2s follow the same Lx - Lmir rela- 
tion is a strong indication that the torus appears optically thin 
at 12.3 //m. 

The correlation we find then requires a dependency of A 
or the covering factor on Lx and thus the bolometric luminos- 
ity LboI- The simplest explanation for this observed effect is 
flaring: If the torus' inner edge is defined by the dust sublima- 
tion radius ro, we obtain A = 27rro • H(ro) where H(r) is the 
scale height of the cloud distribution. For a flaring disk H(r) 
can be written as H(r) oc r" with a being constant. This as- 
sumption yields A oc rl +a . The dust sublimation radius ro fol- 
lows r cc L°.5 oc L° x 5 . This leads to L M ir oc L° x 5+0 - 5a . Our 
observations yield Lmir 00 L l x 6 and thus a x 2. Including NGC 
4579 results in a » 1.5. 

Another possibility is that its scale height depends on 
the radial accretion rate M. The torus model presented by 
iBeckert & Duschll d2004 predicts H{r)Jr oc M 05 . For AGN 
with comparable radiative efficiencies M is proportional to Lx. 
This yields Lmir oc A oc L x 5 which is very close to the correla- 
tion we find. 

The high angular resolution of our MIR observations in- 
dicates that the similarity of the Lmir / Lx ratio for Syls and 
Sy2s is intrinsic rather than caused by extra-nuclear emission. 
Although contributions of dust clumps situated in the ionisa- 
tion cone cannot be ruled out they should not dominate the nu- 
clear MIR flux. Assuming our sample is representative of the 
general AGN population ou r results thus support recent mod- 
els of clumpy dust tori (e.g. iNenkova et ai1l2002l fkonig et alJ 
2006) and, furthermore, indicate that the filling factor of clouds 
is small. The geometry of the torus can be constrained by the 
slope of the Lx - Lmir correlation. Clumpy tori with low filling 
factors might also account fo r the similarity o f Spitzer spectra 
of type 1 and type 2 QSOs dSturm et al.l2006l) . 
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